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ABSTRACT 


Thais Róding, 1798, commonly known as rock-shell, is among the most frequently found gastropod genera 
worldwide on intertidal rocky shores including those of Japan, China, Taiwan and Korea. This group contains 
important species in many marine environmental studies but species-level taxonomy of the group is quite 
complicated due to the morphological variations in shell characters. This study examined the genetic variations 
and relationships among three Korean Thais species based on the partial nucleotide sequences of mitochondrial 
coxl gene fragments. Phylogenetic trees from different analytic methods (maximum parsimony, neighbor-joining, 
and maximum likelihood) showed that 7. bronni and T. luteostoma are closely related, indicating the most 
recent common ancestry. The low sequence divergence found between T. luteostoma and T. bronni, ranging 
from 1.53% to 3.19%, also corroborates this idea. Further molecular survey using different molecular marker 
is required to fully understand a detailed picture of the origin for their low level of interspecific sequence 
divergence. Sequence comparisons among conspecific individuals revealed extensive sequence variations 
within the three species with maximum values of 2.43% in T. clavigera and 1.37% in both T. bronni and T. 
luteostoma. In addition, there is an unexpectedly high level of mitochondrial genotypic diversity within each of 
the three Korean Thais species. The high genetic diversity revealed in Korean Thais species is likely to reflect 
genetic diversity introduced from potential source populations with diverse geographic origins, such as 
Taiwan, Hong Kong, and a variety of different coastal regions in South China and Japan. Additional sequence 
analysis with comprehensive taxon sampling from unstudied potential source populations will be also needed 
to address the origin and key factors for the high level of genetic diversity discovered within the three Korean 


Thais species studied. 
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INTRODUCTION 


Thais Róding, 1798 (commonly known as rock-shell) is one 
of the most frequently found genera of gastropods worldwide 
on intertidal rocky shores including those of Japan, China, 
Taiwan, and Korea. The genus contains about a hundred 
species names (Houart and Gofas, 2010), some of which 
includes predatory species that cause serious damage to the 
oyster culture industry by drilling through the oyster shells 
(Clench, 1947; Brown and Richardson, 1988). The species 
of this genus also play an important role in the marine macro- 
benthic community. They have often been used as an effec- 
tive indicator species in many marine environmental studies 


related to human-induced pollution, such as that from orga- 
notin compounds (Shim et al., 2000; Hung et al., 2001; Tang 
and Wang, 2009) and long-term exposure to heavy metal 
contamination (Rubio et al., 1993; Han et al., 1997; Black- 
more and Wang, 2004). 

Despite their significance in aforementioned environmen- 
tal studies, the species-level taxonomy of this group is quite 
complicated due to morphological variations in shell charac- 
teristics (Tan and Sigurdsson, 1996). Indeed, the morphologi- 
cal features used most frequently for species identification, 
including size, shape and the absence or presence of blotches 
on the nodules of shell surfaces and shell apertures, are high- 
ly variable depending on the local environment (Hayashi, 
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1999; Tan and Liu, 2001). Owing to the substantial variabili- 
ty in shell characteristics, species identification relying on 
shell morphology alone often leads to erroneous results creat- 
ing many synonyms. As a result, the associated nomenclature 
for Thais species is unusually complicated (Clench, 1947; 
Tan, 1995). 

In Korea, Thais species have been included as part of a 
species checklist in the faunal report on mollusks. Earlier 
studies, including some encyclopedias of Korean shells, 
relied entirely on the shell characteristics without careful 
appraisal of species identity (Choe, 1992; Lee and Min, 2002; 
Min, 2004). A taxonomic investigation of Korean muricid 
species based on a morphological comparison of shell and 
radula characteristics showed that three Thais species are 
very abundant in Korean coastal areas (Choe and Park, 1997): 
T. clavigera (Küster, 1858), T. bronni (Dunker, 1860), and 
T. luteostoma (Holten, 1803). In contrast to taxonomic assess- 
ments by most earlier authorities, these species are not very 
morphologically discernible from each other because their 
shell morphology is greatly variable, creating a diverse array 
of local morphotypes according to their geographic origins. 
Therefore, there is considerable taxonomic confusion with 
regard to species-level taxonomy in this group. 

The utility of molecular markers in modern taxonomy has 
been growing rapidly as a reliable tool for examining the 
phylogenetic relationships in different levels of taxonomic 
groups and species identification. As in many other animal 
groups, molecular data from mitochondrial (mt) gene frag- 
ments (e.g., 16S and cox1) provide a wealthy resource for 
assessing the genetic variations and phylogenetic signals 
among closely related species of diverse molluscan taxa 
(Boudry et al., 2003; Park and Kim, 2003; Lam and Morton, 
2006; Reece et al., 2008). The present study examined the 
genetic variations and relationships among three Korean 
Thais species based on partial nucleotide sequences of mt 
coxl gene fragments. In addition, the sequence data presented 
here provides a very useful molecular identification tool for 
the three Korean Thais species, which are often taxonomi- 
cally complicated when based upon the morphological char- 
acteristics alone. 


MATERIALS AND METHODS 


The specimens were collected from the intertidal and/or 
subtidal zones in nine localities of the Korean seashore. 
Voucher specimens used for species identification and subse- 
quent molecular sequencing were deposited in the Marine 
Mollusk Resource Bank of Korea (MMRBK; Chungbuk 
National University, Korea). Identification of the Thais speci- 
mens was carried out based on the shell characteristics, from 


which they were assigned to three species: T. clavigera, T. 
luteostoma, and T. bronni (Fig. 1). A total of 46 specimens 
sampled from nine localities representing three Korean Thais 
species were used for sequencing analysis. The total geno- 
mic DNA was extracted using a tissue kit (Qiagen, Valencia, 
CA, USA) according to the manufacturer's instructions. A 
universal primer set (LCO1490: 5'-GGTCAACAAATCAT 
AAAGATATTGG-3’, HCO2198: 5'-TAAACTTCAGGGT 
GA CCAAAAAATCA-3’) was employed for PCR amplifi- 
cation of the target gene fragment mt cox/ (Folmer et al., 
1994). PCR reactions were performed in a 50 uL reaction 
volume consisting of 10 units of Taq polymerase (Roche, 
Mannheim, Germany), 2.5 mM dNTP mixture, 2.5 mM MgCl; 
and 20 pmole of each primer with the following amplification 
conditions: one cycle of the initial denaturation step at 94^C 
for 2 min, followed by 35 cycles of denaturation at 94°C for 
30 sec, primer annealing at 45°C for 30 sec and elongation 
at 72°C for 1 min with a final extension step at 72°C for 10 
min. The PCR-amplified target gene fragment was purified 
using a QlAquick gel extraction kit (Qiagen) according to 
the manufacturer's protocol. The sequencing reaction was 
performed using a BigDye? terminator v3.1 cycle sequencing 
kit (Applied Biosystems, Foster City, CA, USA), and the 
reaction products were electrophoresed using an ABI 3730XL 
DNA analyzer. Phylogenetic analysis for the mt cox/ dataset 
obtained from the three Thais species was performed using 
different tree-building methods (maximum parsimony [MP], 
neighbor-joining [NJ], and maximum likelihood [ML]) with 
PAUP 4.0b10 (Swofford, 2002) using Thais turbinoides as 
an outgroup. NJ and ML analyses were carried out using the 
HKY +T (0.0913 of gamma shape parameter) model that 
was selected as the ‘best-fit’ substitution model according 
to the Akaike information criterion (AIC) from implementa- 
tion of the Modeltest version 3.7 program (Posada and 
Crandall, 1998). The reliability of grouping in both MP and 
NJ methods was estimated using non-parametric bootstrap 
resampling of 1,000 pseudo-replicates. 


RESULTS AND DISCUSSION 


Phylogenetic relationships among Korean Thais 
species 

The target mt cox] gene fragment determined for all Thais 
species is 658 base pairs (bp) in length with no intra- or 
interspecific length variations. Table 1 lists the GenBank 
accession numbers for each mt genotype and locality data 
for the species. The MP analysis of the mt cox/ dataset show- 
ed two distinct clades with very strong support (Fig. 2): one 
representing the assemblage of mt genotypes from T. clavi- 
gera and the other containing two groups of mt genotypes, 
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Fig. 1. Shell morphology of three Korean Thais species. T. clavigera (1a, 1b), T. bronni (2a-2f), T. luteostoma (3a-3d). Scale bars=2 cm. 
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Table 1. Geographic origin, sampling locality and GenBank accession numbers for Thais species used in this study 


No. of individuals GenBank 


Species Geographic origin examined Genotype accession no: 

Thais luteostoma Tonggumi, Ulleung-gun 10 TLO1 (1) HQ852742 
Gyeongsangbuk-do TLO2 (1) HQ852743 

TLO3 (1) HQ852744 

TL04 (3) HQ852745 

TLO5 (1) HQ852746 

TLO6 (1) HQ852747 

TL10 (1) HQ852751 

TL12 (1) HQ852753 

Daeyeong, Jeju-do 6 TLO2 (1) HQ852743 

TLO4 (1) HQ852745 

TLO7 (1) HQ852748 

TLO8 (1) HQ852749 

TLOO9 (1) HQ852750 

TL11 (1) HQ852752 

Seogwipo-si, Jeju-do 2 TL10 (2) HQ852751 

Thais bronni Seogwipo-si, Jeju-do 4 TBO1 (1) HQ852754 
TBO2 (1) HQ852755 

1803 (1) HQ852756 

1804 (1) HQ852757 

Dodong, Yokji-myeon, 9 1803 (1) HQ852756 

Tongyeong-si, 7805 (1) HQ852758 

Gyeongsangnam-do TBO6 (1) HQ852759 

TBO7 (1) HQ852760 

TBO8 (1) HQ852761 

TBO9 (1) HQ852762 

TB10 (2) HQ852763 

TB11 (1) HQ852764 

Gajin-ri, Goseong-gun, Gangwon-do 2 TBO7 (1) HQ852760 

TBO8 (1) HQ852761 

Wonbuk-myeon, Taean-gun, 2 TB12 (1) HQ852765 

Chungcheongnam-do TB13 (1) HQ852766 

Thais clavigera Ganwoldo-ri, Seosan-si, Chungcheongnam-do 1 TCO1 (1) HQ852767 
Hoenggando, Chujado, Jeju-do 1 TCO2 (1) HQ852768 

Dodong, Yokji-myeon, 2 TC03 (1) HQ852769 

Tongyeong-si, Gyeongsangnam-do 7004 (1) HQ852770 

Geoado-ri, Taean-gun, 7 TC03 (1) HQ852769 

Chungcheongnam-do 70604 (1) HQ852770 

70605 (1) HQ852771 

70606 (1) HQ852772 

1Շ07 (1) HQ852773 

TCOS8 (1) HQ852774 

7009 (1) HQ852775 

Gajin-ri, Goseong-gun, Gangwon-do 1 TC10 (1) HQ852776 

Thais haemastoma 1 EU073051 
Thais turbinoides South of New Caledonia 1 HQ852777 


representing T. luteostoma and T. bronni. The grouping pat- 
terns among the three species in the MP analysis were iden- 
tical to those of the NJ and ML analyses (not shown). The 
bootstrap values for these two clades were very high (10096 
for T. clavigera in both MP and NJ methods, and 96% and 
10096 for T. bronni-T. luteostoma in MP and NJ analyses, 
respectively). The branches representing 7. luteostoma (9296 
and 96% bootstrap supports in MP and NJ, respectively) and 
T. bronni (88% and 81% supports in MP and NJ analyses, 


respectively) also received relatively strong support. In con- 
trast to the well-structured relationships, the inferred tree 
showed no particular grouping pattern among the mitochon- 
drial genotypes within each of the three species. Irrespective 
of the analytical methods (MP, NJ, and ML), the reconstruct- 
ed phylogenetic trees consistently depicted 7. bronni and T. 
luteostoma as being more closely related to each other than 
to T. clavigera. This sister group relationship was robustly 
supported by a bootstrap value of 9696 (MP) and 10096 (NJ) 
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Fig. 2. Bootstrapped consensus tree of 28 equally parsimony trees from heuristic analysis for cox1 sequence data of three Korean 
Thais species. Bootstrap supporting values from the maximum parsimony method are indicated above the internal branches if >50%. 
The bootstrap values for each node from neighbor-joining analysis using the HKY+I (0.0913 of gamma shape parameter) model as 
the ‘best-fit’ substitution model selected under the Akaike information criterion from the implementation of Modeltest 3.7 program 


(Posada and Crandall, 1998) were also shown below the branches. 


in the different analytic methods (MP and NJ). This result 
corresponds to the morphological similarities revealed by 
some earlier studies using comparative analysis of morpho- 
logical characters (shell, radula, and penis morphology etc.) 
(Choe and Park, 1997; Tan and Liu, 2001). 

Aside from the shell external characteristics, traditional 
taxonomy has appreciated the radula morphology as one of 
the most useful features for the classification of a wide range 
of gastropod groups, including muricid species (Cooke, 1919; 
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Wu, 1965). Nevertheless, the radula characteristics are often 
quite similar among closely related congeneric species, 
making species identification rather difficult based on the 
radular traits alone (Fujioka, 1985; Tan and Liu, 2001). Some 
earlier authorities reported that the rachidian teeth of T. 
luteostoma and T. bronni were quite similar but noticeably 
different from those of T. clavigera in both the shape and 
number of denticles. There are usually more than six denticles 
between the lateral and marginal cusps recorded for the 
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Korean population of T. luteostoma and T. bronni (6 and 7- 
8, respectively) (Choe and Park, 1997). Tan and Liu (2001) 
observed some variations from Hong Kong population of T. 
luteostoma in the denticle number, ranging from five to eight. 
The rachidian tooth of T. clavigera has fewer denticles (nor- 
mally 3-4 denticles between lateral and marginal cusps). 
Although variations in radula characteristics according to 
the local populations are not uncommon, particularly in the 
denticle number, similarities in the radula characteristics 
between T. bronni and T. luteostoma are in line with the 
close relationship uncovered by genetic analysis of the mol- 
ecular data in this study. The closer relationship between 7. 
bronni and T. luteostoma is also supported by a comparison 
of the penis morphology. Indeed, the morphology of the 
male reproductive organ is generally considered species- 
specific and appears to be very useful for species identifica- 
tion in this genus. Exceptions were found for T. luteostoma 
and T. jubilaea, whose penis morphology was more or less 
similar (Tan and Sigurdsson, 1990; Tan and Liu, 2001). Inter- 
estingly, the penis morphology of 7. luteostoma was quite 
different from that of T. clavigera, but similar to that of T. 
bronni from the Korean population. T. bronni has a recurved 
hooklet-shaped penis with a gradually tapered terminal end 
(Park, 1996) and these characteristics were also found in 
Hong Kong (Tan and Liu, 2001) and Korean (personal obser- 
vations) populations of T. luteostoma. The morphological 
similarities of both radula and penis are consistent with the 
current molecular data, indicating that T. luteostoma and T. 
bronni are very closely related. The inclusion of broader 
taxon sampling will be needed to place the Korean population 
of species within the phylogenetic framework of the genus. 


Mt DNA variations in Korean Thais species 

A total of 35 mt genotypes were discovered from 46 indivi- 
duals of nine local populations of three Korean Thais species 
(Aligned cox] sequences for Thais species are shown in Su- 
pplementary data, Appendix 1). Each of these three species 
included a large number of mt genotypes: 10 for T. clavigera, 
13 for T. bronni and 12 for T. luteostoma. Table 2 lists the 
uncorrected pairwise (p) distance among mt genotypes. The 
sequence divergence of T. turbinoides (outgroup species) 
from Korean species ranged from 15.2% (between TLO8 of 
T. luteostoma and T. turbinoides pair and between TBO7 of 
T. bronni and T. turbinoides pair) to 16.9% (between TCO1 
of T. clavigera and T. turbinoides pair). The highest interspe- 
cific sequence divergence (17.796) was observed between T. 
haemastoma and the genotype TL11 from T. luteostoma. On 
the other hand, sequence comparison within Korean congener- 
ic species revealed the 7. clavigera sequences to differ con- 
siderably from those of T. luteostoma and T. bronni, with 
maximum differences ranging from 7.14% (between TCOI 


and TLO7 pair) to 7.75% (between TCO1 and TBO6 pair). In 
contrast, interspecific sequence divergence between T. luteo- 
stoma and T. bronni was relatively very low, ranging from 
1.52% (for TL10 and TBO3 pair) to 3.19% (for TLO5 and 
TBO6 pair). The lowest interspecific difference value (1.52%) 
detected between T. luteostoma and T. bronni was even 
lower than the highest divergence value detected within T. 
clavigera (2.43% between TCO1 and TCO08). This low diver- 
gence value found between T. luteostoma and T. bronni 
implies that cladonegenetic split of these two species is very 
young with a very recent common origin. Further molecular 
survey using different molecular marker is required to fully 
understand a detailed picture of the origin for their low level 
of sequence divergence. In addition, sequence comparison 
of conspecific individuals revealed extensive sequence varia- 
tions within each of the three species. The intraspecific differ- 
ences among the 7. clavigera mt genotypes were remarkably 
high, with a maximum 2.43% divergence of a pairwise com- 
parison (found between TCO1 and TC08). The sequence 
divergence within species in each of T. bronni and T. luteo- 
stoma was also high: the maximum value was detected bet- 
ween TB06-TB02 and TB06-TBO7 pairs in T. bronni (1.3796) 
and between TL05-TL06, TLO5-TLO7, and TLO05-TL 12 pairs 
in T. luteostoma (1.37%). 

Along with the considerable genetic divergence within 
species, there were a large number of mt genotypes discover- 
ed in each of the three species. Thirty five mt genotypes were 
detected from 46 individuals of the three species sequenced. 
Of these, a small fraction of mt genotypes (9 of 35; 25.796 
of total genotypes detected) is represented by more than one 
individual (Table 1, Fig. 2). With the exception of these, all 
other mt genotypes (26 of 35; 74.396 of total genotypes de- 
tected) were presented by a single individual. These exclu- 
sive mt genotypes were shown to be unresolved polytomies 
at the terminal tips in the reconstructed phylogenetic tree 
(Fig. 2). This unexpectedly high genotypic diversity encoun- 
tered within each of the three Korean Thais species is assum- 
ed to be a reflection of the genetic diversity introduced from 
potential source populations with diverse geographic origins, 
such as Taiwan, Hong Kong, and many different coastal 
regions of South China and Japan. Although a large expansion 
of deep oceanic water is considered an obstacle to the long- 
distance dispersal of many marine benthic taxa, some coastal 
invertebrate taxa have crossed this barrier using a prolonged 
pelagic larval form (Scheltema, 1971, 1986, 1988) or by 
rafting (passive transportation of sessile-form juveniles atta- 
ched to drifting objects) (ÓFoighil et al., 1999). In cases of 
long range dispersion, the juvenile forms of these animal 
groups may have the potential to transport and establish 
successful colonies thousands of kilometers away in down- 
stream habitats that are geographically distant from up- 
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stream source populations using the oceanic current sys- 
tems. Indeed, Thais species have a wide distribution range 
in the Asian Pacific Ocean from the South China Sea area 
(including Hong Kong, Taiwan) to Northeast Asia (includ- 
ing Japan and Korea) (Kuroda et al., 1971; Choe and Park, 
1997; Tan, 2000; Tan and Liu, 2001). It is possible to envis- 
age that a continuous influx of larval immigrants from their 
potential source populations might contribute to some ex- 
tent to the near-shore malacofauna of the Korean Peninsula. 
Of the western boundary currents in the North Pacific 
Ocean, the Kuroshio Current System is a major component 
of the most influential warm current that leaves the east coast 
of Taiwan. Its surface water effects travel northeastward to 
many regional rocky shores in eastern China, mainland Japan, 
and the Korean Peninsula by streams of the Tsushima Cur- 
rent and Yellow Sea Current. These current systems running 
northeastward may be able to transport planktotrophic larvae 
or drifting objects that carry sessile forms of life (e.g., egg 
capsules of Thais species) to many different downstream 
islands and/or continental rocky shores. The routine influx 
of these types of immigrants from geographically distant 
source populations with different geographic origins (e.g., 
Taiwan, Hong Kong, and many different coastal regions of 
South China and Japan) and their ability to colonize in down- 
stream habitats might act as an effective supplier of new 
recruitment, thereby shaping the contemporary population 
structure of Korean Thais species. Under this assumption, the 
mt genotypic diversity found among Korean Thais popula- 
tions is likely to reflect a subset of the genetic diversity of 
their potential source populations. Sampling from their po- 
tential source populations was not included in the present 
study. Therefore, further comprehensive taxon sampling from 
unstudied potential source populations will be needed to de- 
termine the origin and key factors responsible for the high 
level of genetic diversity discovered within each of the three 
Korean Thais species examined. 
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